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A method of processing seismic data 



The present invention relates to a method of processing seismic data and, in particular, 
to a method of processing marine seismic data to reduce the effects of "ghost 
reflections". The present invention can be applied to processing both marine seismic 
data acquired in a flat sea and marine seismic data acquired in a rough sea. The present 
invention also relates to processing seismic data so as to estimate the depth of a seismic 
receiver, or so as to monitor the accuracy of a seismic receiver. 

Figure 1 is a schematic diagram of a marine seismic survey in which seismic energy is 
emitted from a source 1 and detected by a seismic receiver 2 at a depth h below the 
surface 6 of the sea. Energy emitted from the source is reflected by the seabed 3 or by a 
reflector 4 below the seabed 3 and is then detected by the receiver. This path of seismic 
energy is labelled 5 in Figure 1. Information about the geological structure of the 
earth's interior can be derived from the reflected seismic energy incident on the receiver 

The seismic receiver 2 shown in Figure 1 is a streamer, which is a type of seismic 
receiver that is frequently used in marine seismic surveying. A streamer contains a 
plurality of sensors Si, S2,...Sn such as pressure sensors and/or particle velocity sensors 
. distributed along its length, which may be some hundreds of metres, and is thus able to 
measure the reflected seismic energy at a number of points simultaneously. A streamer 
is generally suspended from one or more floats 8, so that all the receivers of the 
streamer are at the same depth in a flat sea. 

In addition to the desired path 5 of seismic energy shown in Figure 1, other seismic 
energy paths will occur as a result of seismic energy being reflected or scattered from 
the sea surface 6. These paths are known as "ghost reflections". For example, reference 
7 in Figure 1 shows a ghost reflection in which the seismic energy reflected by the 
reflector 4 is not directly incident on the receiver 2, but undergoes a further reflection at 
the sea surface 6 before reaching the receiver. Down-going sea-surface ghost 



WO 02/01254 



PCT7IB01/01150 



2 

reflections are an undesirable source of contamination of seismic data, since they 
obscure the interpretation of the desired up-going reflections from the earth's interior. 

Ghost reflections also arise from sea surface reflections or scattering near the source 1. 
The overall seismic energy signal received at a point below the source will contain a 
component that has come direct from the source and a component that has come via a 
reflection or scattering at the sea surface. 

Removing the ghost reflections from seismic data is for many experimental 
configurations equivalent to up/down wavefield separation of the recorded data. In such 
configurations the down-going part of the wavefield represents the ghost signal and the 
up-going wavefield represents the desired signal. The up/down wavefield separation 
can be achieved by applying a suitable filter to the seismic data. 

The present invention provides a method of de-ghosting seismic data, the method 
comprising the steps of: 

a) obtaining seismic data using a seismic source and a seismic receiver disposed 
within a fluid medium; 

b) determining a derivative with respect to time of the fluid pressure at a point in 
the vicinity of a selected one of the seismic source and the receiver; 

c) detennining a derivative with respect to a first direction of the fluid pressure at a 
point in the vicinity of the selected one of the seismic source and receiver; 

d) determining a derivative with respect to a second direction of the fluid pressure 
at a point in the vicinity of the selected one of the seismic source and receiver from the 
results of steps (b) and (c); and 

e) processing the seismic data using the derivative with respect to a second 
direction of the fluid pressure, thereby to reduce the effects in the processed seismic 
data of seismic energy reflected and/or scattered at the sea surface. 

In a typical marine seismic survey, the time-derivatives of the pressure and the 
derivatives in one direction, for example the horizontal direction, are readily available, 
whereas the derivatives in other directions, for example the vertical direction, are not 
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directly known and have to be estimated. The present invention uses the readily 
available derivatives to provide a better estimate of a derivative that is not directly 
known, and this allows more accurate de-ghosting of the seismic data. 

In a preferred embodiment step (a) comprises obtaining seismic data at two different 
depths below the surface of the fluid medium, and in a particularly preferred 
embodiment step (a) comprises measuring the fluid pressure at two different depths 
below the surface of the fluid medium. This allows the derivative of the de-ghosting 
step to be further improved. 

In an alternative preferred embodiment, step (a) comprises measuring the fluid pressure 
and the particle velocity. If the depth is known, this embodiment allows the pressure 
and particle velocity sensors to be monitored. For example, the particle velocity can be 
computed from the depth and the measured pressure values, and the computed value 
may then be compared with the measured values of the particle velocity. Alternatively, 
the depth may be estimated from the measured values of the particle velocity and 
pressure. 

A second aspect of the present invention provides a method of calibrating a seismic 
receiver as claimed in claim 20. This method enables the accuracy of measurement of 
particle velocity to be monitored. 

A third aspect of the present invention provides a method of calibrating a seismic 
receiver as claimed in claim 22, and a fourth aspect of the invention relates to a method 
of determining the depth of a seismic receiver as claimed in claim 24. These aspects of 
the invention relate to a receiver that can measure both the particle velocity and the fluid 
pressure. 

Other features of the invention are set out in the dependent claims. 

Preferred embodiments of the present invention will now be described by way of 
illustrative example with reference to the accompanying figures in which: 
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Figure 1 is a schematic view of a marine seismic survey illustrating the production of 
ghost reflections; 

Figure 2 is a schematic profile view of a rough sea surface; 

Figure 3(a) shows the results of a prior art de-ghosting method for a streamer depth of 
3.5m in a flat sea; 

Figure 3(b) shows the results of a prior art de-ghosting method for a streamer depth of 
5.5m in a flat sea; 

Figure 4(a) shows the results of a de-ghosting method according to a first embodiment 
of the present invention for a streamer depth of 3.5m in a flat sea; 

Figure 4(b) shows the results of a de-ghosting method according to the first embodiment 
of the present invention for a streamer depth of 5.5m in a flat sea; 

Figure 5(a) shows the results of a prior art de-ghosting method for a streamer depth of 
3.5m in a rough sea; 

Figure 5(b) shows the results of a prior art de-ghosting method for a streamer depth of 
5.5m in a rough sea; 

Figure 6(a) shows the results of a de-ghosting method according to the first embodiment 
of the present invention for a streamer depth of 3.5m in a rough sea; 

Figure 6(b) shows the results of a de-ghosting method according to the first embodiment 
of the present invention for a streamer depth of 5.5m in a rough sea; 

Figure 7(a) shows the results of a de-ghosting method according to a second 
embodiment of the present invention for a streamer depth of 3.5m in a flat sea; and 
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Figure 7(b) shows the results of a de-ghosting method according to the second 
embodiment of the present invention for a streamer depth of 3.5m in a rough sea. 

The equation relating pressure p to particle velocity v in an acoustic medium with a 
bulk modulus K is: 



P = -KV.v (1) 



In equation (1), p denotes the time-derivative of p. The pressure p and particle 
velocity v are also related by the equation of motion, which has the following form 



v=--V/> (2) 
P 



Here, pis the density of the acoustic medium. The equation of motion (2) and equation 
(1) together complete the set of first order partial differential equations describing 
acoustic wave propagation. These equations can be combined to give the well-known, 
second-order, acoustic wave equation: 



p = a 2 V 2 p (3) 



In equation (3), a = ^Kl p is the wave propagation speed in the acoustic medium. 



Since the curl of the gradient of an acoustic wavefield is zero (this being a fundamental 
vector relation), equation (2) yields: 



V x v = 0 (4) 



A preferred embodiment of the invention relating to de-ghosting a streamer located in 
the vicinity of the sea-surface will now be described. In this embodiment the streamer 



WO 02/01254 



PCT/IB01/01150 



6 

is provided with a plurality of pressure sensors distributed along its length. Although 
the propagation of seismic energy is three-dimensional in any real seismic surveying 
arrangement, this example will consider only 2D wave propagation in the x-z plane (z 
increases downwards, as indicated in Figure 1) for clarity of explanation. The invention 
can, however, be applied to full three-dimensional propagation of seismic energy. 

It has been reported, for example in co-pending UK patent application No 9906456.0, 
that de-ghosting of seismic data can be carried out if the vertical derivative of pressure 
d z p is known. Thus, it is necessary to compute the vertical derivative of pressure 

along the streamer at a point (x 0 , zq). This can be achieved by writing an estimate of the 
vertical derivative of pressure as a Taylor expansion in the ^-direction of p around (x# g 
zq). After rearranging the terms of the Taylor expansion we obtain: 

a,^o^o)- ^ 0>Zo) "f (xo>2 °" fe) + f^o^o) + ^) (5) 

In equation (5), h is the vertical distance of the streamer at point (xo,zo) below the sea- 
surface, d 2 denotes a first-order derivative in the z-direction and d a denotes a second 

order derivative. The OQi ) term shows that the largest error terms have factor h in 
them and can therefore be ignored to a reasonably accurate approximation. 

In the case of a marine seismic survey that has a single streamer it is not possible to 
measure d a p directly. Pressure is measured only at the streamer, and this may be at a 
constant depth in a flat sea. However, the second-order derivative of/? in the vertical 
direction, d a p(x 0 9 z 0 ) 9 can be related to the time derivative of pressure and to the 

derivative of pressure in a direction along the streamer, which is assumed to be 
horizontal, by using the wave equation. This allows us to increase the accuracy of a 
derivative estimate. A correction of this type is known as a Lax-Wendroff correction in 
the synthetic finite-difference literature, and is described by Lax and Wendroff in 
'Difference Schemes for Hyperbolic Equations with high order accuracy", Comm. Pure 
Appl. Math Vol. 27 (1964), and by Blanch and Robertsson in "Geophysical Journal 
International", Vol. 131,pp381-386 (1997). 
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Thus, according to the invention Equation (3) is used to convert derivatives that are not 
known into derivatives that can be computed from the measured pressure data. From 
equation (3) the following relationship can be obtained: 

9 a p(x 0 9 z 0 ) = ~ p(x 0 , z 0 ) - d , z 0 ) (6) 

a 

By substituting the above formula for d^p, Equation (5) then becomes: 

^ 0 ,zJ = ^^ + ^^ OJ zJ-|5„(x 0 ,zJ + C?(AO (7) 

Here we have made use of the fact that the pressure vanishes at the sea-surface: 
pfao, zo-h) = 0. 

The significance of equation (7) is it does not contain d a p. In equation (7), d z p is 

expressed in terms of the derivative of the pressure with respect to time and the 
derivative of the pressure with regard to the ^-direction (the horizontal direction). Both 
of these quantities can be readily estimated from the available data. The time-derivative 
of the pressure can be readily computed from pressure measurements made by a sensor 
over time. Moreover, as noted above a streamer contains a plurality of pressure sensors 

* 

disposed along its length, and the derivative of the pressure in the ^-direction can be 
readily computed from pressure measurements made at a given time by adjacent 
pressure sensors along the streamer. 

Equation (7) also provides a means of computing both the x-component and the z- 
component of particle velocity along the streamer, using equation (2), provided that the 
streamer depth at each point along the streamer is known at all times. 

The derivatives with respect to x, d x p , in the expression for derived from equation 
(2) and d^p in equation (7), can be computed by differentiating pressure measurements 
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along the streamer. As described in co-pending UK patent application No 9906456.0, 
computing d x p enables the seismic data to be de-ghosted. Applying the following 



The variable p denotes the recorded up-going pressure (deghosted), i is the square root 



of — 1, k z = y(W a) 2 - k\ is the vertical wavenumber, is the horizontal wavenumber, 
and w is the angular frequency. 

The method for estimating d z p according to the invention in principle is independent 
of the roughness of the sea, as is the deghosting step described in equation (8). 

Thus, according to the present invention, the derivative in the z (vertical) direction of 
the pressure is determined from the pressure measurements made along a horizontal 
line. Once the derivative of the pressure in the z direction has been calculated, the 
seismic data can be de-ghosted using equation (8). 

Figures 3(a) to 6(b) compare results of the present invention with results obtained by a 
prior art method. These figures are simulations obtained using a flat sea surface or 
using a rough sea-surface having the profile shown in Figure 2. It will be seen that the 
rough sea surface of Figure 2 has a significant wave height of around 2m. In the 
simulation it was assumed that the pressure is being recorded by a horizontal streamer at 
the points marked with an asterix, namely at x = -12.5m, x = -137.5m, x = -262.5m and 
x - -387.5m. The source is a 50Hz zero-phase explosive source located at x = 0 and at a 
depth of 50m. 

Figures 3(a) and 3(b) show a comparison between d 2 p calculated explicitly using a 

finite - difference method (solid line) and an estimate calculated without a Lax- 
Wendroff correction (broken lines). The traces in each figure correspond to recording 



filter to the seismic data provides de-ghosting, valid for all incidence angles: 




(8) 
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locations x = -12.5m (uppermost trace), x = -137.5m, x = -262.5m, and x = -387.5m 
(lowest trace). 

In Figures 3(a) and 3(b) d z p was calculated simply by dividing the pressure recorded 

by the sensor on the streamer by the vertical distance between the recording location 
and the sea-surface (as noted above, the pressure becomes zero at the sea-surface). 
Thus, this approximation assumes that d z p is constant between the sea-surface and the 

streamer. Figure 3(a) shows the results of a streamer depth of 3.5m, and Figure 3(b) 
shows the results for a streamer depth of 5.5m. It will be seen that computing d z p by 

the approximation of dividing the recorded pressure at the streamer by the vertical 
distance between the recording location and the sea-surface is not particularly accurate. 

Figures 4(a) and 4(b) again relate to a flat sea and show (in dashed lines) results of a 
method of the present invention in which d z p was calculated using equation (7) above 

compared with the results (solid lines) of computing d z p using a finite-difference 

calculation. It will be seen that this method of the present invention provides a much 
more accurate estimate of d 2 p than does the simple geometric averaging used in 

Figures 3(a) and 3(b). 

Figures 5(a) to 6(b) show results of computing d z p for the rough-sea shown in Figure 

2. Figures 5(a) and 5(b) show results for the prior art method. In these figures, the 
results in broken lines are for a calculation of d z p by the prior art method of dividing 

the pressure recorded by the sensor on the streamer by the vertical distance between the 
recording location and the sea-surface, while the results in solid lines are a finite- 
difference calculation for d z p . 

Figures 6(a) and 6(b) show results obtained by using the present invention. In these 
Figures, the solid lines again show the result of a finite-difference calculation for d z p , 
whereas the broken lines show the results of computing d z p using equation (7). It will 
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again be seen that the method of the present invention provides much more accurate 
estimates of d z p than does the simple geometric averaging method. 

In each of Figures 4(a), 4(b), 5(a), 5(b), 6(a) and 6(b), the traces correspond to recording 
locations x = -12.5m (uppermost trace), x = -137.5m, x = -262.5m, and x = -387.5m 
(lowest trace). 

Since the present invention provides more accurate values of d z p , the accuracy of the 
de-ghosting using the filter of equation (8) will also be improved. 

In a further embodiment of the invention, the accuracy of the estimates of d 2 p is further 

increased by eliminating the next higher order term in the Taylor expansion, and 
increasing the accuracy to 0(h 3 ). This requires solving an implicit system of equations 
for d z p at each frequency: 



One way of solving equation (9) is in the f-x domain using a second order accurate 
derivative estimate to approximate the term on the left-hand side. It may also be 

preferable to find filter approximations to equation (9) since equation (9) has effectively 
introduced a pole that could cause instabilities or amplify noise. 

Figures 7(a) and 7(b) shows the results of this higher-order de-ghosting method of the 
present invention. These figures show the result of calculating d z p using a finite- 
difference method (solid lines) and the results according to the invention of calculating 
d z p using equation (9) (broken lines). The four traces in each figure again correspond 

to recording locations of x = -12.5m (uppermost trace), x = -137.5m, x = -262.5m, and x 
= -387.5m (lowest trace). Figures 7(a) and 7(b) both relate to a streamer depth of 3.5m, 
and Figure 7(a) shows results for a flat sea and Figure 7(b) shows the results for the 



d x 1- 



r 




-f «»J<*.*o) + 0(h>) (9) 
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rough sea surface of Figure 2. It will be seen that estimating d z p using equation (9) 
produces generally good results. 

The present invention is applicable to a seismic receiver in the vicinity of the sea- 
surface, such as a streamer that is towed at a shallow depth below the sea-surface. The 
maximum allowable depth of the streamer below the sea-surface is generally similar to 
the maximum separation used to separate the two streamers in a twin streamer 
configuration, since the largest inaccurate term in both cases is of 0(h ) 9 or of 0(h ) if 
equation (9) is used. It is preferable that the depth of the streamer below the sea surface 
does not exceed half the wavelength of seismic energy emitted by the source. 

The present invention requires that the vertical distance below the sea surface of each of 
the sensors distributed along the streamer must be known. A conventional streamer is 
generally provided with a depth sensor so that an operator can monitor the depth of the 
sensor. Provided that this depth sensor is sufficiently accurate, it will be possible to use 
the average depth of the streamer as measured by this sensor as the depth of each of the 
recording elements. This should provide good results in flat sea conditions. For 
improved accuracy in rough sea conditions, however, it is preferable that a plurality of 
depth sensors, such as hydrostatic sensors, are provided on the streamer, so that the 
distance of each pressure sensor below the instantaneous sea-surface can be recorded. 
Alternatively, it is possible that the depth of each sensor can be estimated from the 
acquired seismic data. 

The embodiments of the present invention described above relate to a seismic receiver 
that is provided with sensors for measuring the pressure p of incident seismic energy. 
As previously noted, however, it is also possible to provide a seismic receiver with 
sensors for measuring particle velocity v. An alternative embodiment of the invention 
is suitable for use with a seismic receiver, for example a streamer, that has both pressure 
sensors and particle velocity sensors. Such a streamer can directly measure both the 
fluid pressure p and the particle velocity v. The quantities p and v are not independent 
from one another, however, but are related by equation (1) or (2), so that there is an 
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extra degree of freedom in the subsequent processing of the seismic data acquired by the 
streamer. 

In one embodiment of the invention, it is assumed that the depth of the streamer is 
known. This embodiment allows the output from the particle velocity sensors to be 
compared with the output from the fluid pressure sensors. This comparison can be 
made, for example, by calculating the particle velocity from the measured values for 
pressure and from the depth of the streamer, and comparing these calculated values of 
the particle velocity with the directly measured values for the particle velocity. These 
values should, of course, be the same if the particle velocity and the fluid pressure have 
been measured accurately. Alternatively, the pressure can be calculated from the 
measured values for the particle velocity and from the depth of the streamer, and these 
calculated values of the pressure can be compared with the directly measured values for 
the pressure. This process can be considered as calibrating the pressure sensors relative 
to the particle velocity sensors, and allows the accuracy of the sensors to be monitored. 

This process for monitoring the accuracy of the sensors can be carried out while 
acquired seismic data is being processed by a de-ghosting method of the present 
invention. Alternatively, it can be carried out on its own, for example to check the 
sensors at the start of a seismic survey. This allows the extra degree of freedom 
provided by measuring both p and v to be used to calibrate the sensors. 

In this embodiment, the depth of the streamer can be estimated, for example from 
knowledge of the arrangement for suspending the streamer from floats. For improved 
accuracy, however, it is preferable that the depth of the receiver is measured directly. 

In another embodiment of the invention, the extra degree of freedom resulting from the 
measurement of both p and v is used to calculate the distance between each pressure 
sensor and the sea-surface. In the case of a flat sea, this gives the depth and shape of the 
streamer. 
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A further embodiment of the invention again relates to monitoring the accuracy of a 
streamer that is provided with both fluid pressure and particle velocity sensors. In the 
embodiment the accuracy of the particle velocity sensors is monitored by using equation 
(4), which states that the curl of the particle velocity should be identically zero. This in 
principle allows the particle velocity sensors to be monitored since, if v has been 
measured accurately, then calculating Vxv should give zero so that Vxv acts as a 
measure of the quality of the measurement of v. In practice, however, it has hitherto 
been difficult to compute Vxv with sufficient accuracy for this to be a reliable method 
of monitoring the particle velocity sensors. Measurements of v are available at many 
points along the streamer but are only available at one depth, and this makes it difficult 
to compute the derivatives with respect to depth of components of v . 

According to this embodiment of the invention, a Lax-Wendroff correction is used to 
improve the accuracy of the computation of the derivatives with respect to depth of the 
components of v , and thereby improve the accuracy with which Vxv can be 
computed. The wave equation (3) can alternatively be written in terms of v rather than 
p 9 and using a Lax-Wendroff correction enables the derivative with respect to depth of 
the components of v to be expressed in terms of derivatives with respect to time and 
derivatives with respect to distance along the streamer. This enables the derivative with 
respect to depth of the components of v to be computed with greater accuracy; this is 
analogous to the increase in accuracy of computing d a p in previous embodiments. 

A further embodiment of the invention is for use with a configuration in which two 
streamers are towed one above the other in the vicinity of the sea-surface. In this 
configuration, the pressure is measured at the level of each streamer, and it is also 
known that the pressure at the sea-surface is zero. This receiver configuration therefore 
provides sufficient information to enable the second-derivatives of pressure in both the 
horizontal direction and the vertical direction to be calculated from the measured 
pressure data. It is thus possible to determine d z p directly from equation (5), since all 

the terms on the right-hand side of equation (5) can be computed directly from the 
measured pressure data. In this embodiment, the wave equation (3) which relates 
spatial derivatives of/? to time derivatives of/?, is not required in the computation of 
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d z p , and can be used to provide a constraint on the streamer depth or on the streamer 

separation in a de-ghosting step, since limits on the streamer depth or streamer 
separation can be determined. 

Alternatively in this embodiment of the invention, a Lax- WendrofF relation can be used 
to increase the accuracy of the estimate of the vertical derivative between the streamers. 
This will allow the streamers to be towed further apart from each other while properly 
handling higher frequencies in the data. For this situation, the depth derivative is 
centred at a depth midway between the depths of the two streamers so that the Taylor 
expansion of the approximation takes the form: 

p(x '- H + hn \^'-*° - hn) -^a^„z.) + oiH') (io) 

In Equation (10), h is the vertical separation between the two streamers. 

Again, by using equation (3) to write d^p in terms of derivatives with respect to time 

or the x-direction, it is possible to obtain an implicit system of equations to solve for 
d z p. 

The embodiments of the present invention described above have related to the de- 
ghosting the signal recorded at the receiver. The invention can, however, also be used 
for source-side de-ghosting. In principle this could be done by directly measuring the 
fluid pressure at the source, but in practice it is possible to make use of the principle of 
reciprocity to avoid the need to do this. 

The principle of reciprocity is a fundamental principle of wave propagation, and states 
that a signal is unaffected by interchanging the location and character of the sources and 
receivers. For example, if a surveying arrangement with a seismic source at point A and 
a receiver array at point B gives a certain signal at the receiver array, then using a single 
receiver at point A and a source array at point B would lead to the same signal, provided 
that the source array corresponds to the receiver array. (By "corresponds", it is meant 
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that the source array contains the same number of sources as the receiver array has 
receivers, and that the sources in the source array are arranged in the same positions 
relative to one another as the receivers in the receiver array). 

In consequence of the principle of reciprocity, the signal emitted by an array of seismic 
sources can be de-ghosted by making measurements using a receiver array, whether the 
seismic sources are actuated simultaneously or consecutively. The recorded signal is 
analogous to that recorded in a reciprocal arrangement where the source and receiver 
locations are interchanged. The method outlined hereinabove can therefore also be 
applied to the source array that it is desired to de-ghost. The signal produced at the 
receiver array by an array of a plurality of seismic sources is measured, and a de- 
ghosting filter is derived from this measured signal. By the principle of reciprocity, this 
filter can be used to de-ghost the signal emitted by the source array. 

The present invention can provide good results in source-side de-ghosting, since the 
depth of the source(s) in a typical marine seismic survey is generally well-known, for 
example because the source(s) is/are suspended from a boat by a known length of cable. 
The present invention can particularly be applied to a full three-dimensional VSP 
(vertical seismic profiling) acquisition pattern, where a large number of seismic sources 
are generally employed in a relatively small area. The full three-dimensional version of 
the present invention can then be applied efficiently. 

The methodology of the present invention can also be used to design seismic source 
arrays that have two or more marine vibrators, air guns, or water guns. 
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CLAIMS: 



1. A method of de-ghosting seismic data, the method comprising the steps of: 

a) obtaining seismic data using a seismic source and a seismic receiver disposed 
within a fluid medium; 

b) determining a derivative with respect to time of the fluid pressure at a point in 
the vicinity of a selected one of the seismic source and the receiver; 

c) determining a derivative with respect to a first direction of the fluid pressure at a 
point in the vicinity of the selected one of the seismic source and receiver; 

d) determining a derivative with respect to a second direction of the fluid pressure 
at a point in the vicinity of the selected one of the seismic source and receiver from the 
results of steps (b) and (c); and 

e) processing the seismic data using the derivative with respect to a second 
direction of the fluid pressure, thereby to reduce the effects in the processed seismic 
data of seismic energy reflected and/or scattered at the sea surface. 

2. A method as claimed in claim 1 wherein the derivative of the fluid pressure with 
respect to the second direction is determined from the derivative with respect to time of 
the fluid pressure and from the derivative with respect to the first direction of the fluid 
pressure using the equation 

P = & V p 

where p is the fluid pressure, a = ^Kl p where K is the bulk modulus and p is the 
density of the fluid, and p denotes the second-order time-derivative of p. 

3. A method as claimed in claim 1 or 2 wherein step (d) comprises determining a 
first-order derivative of the pressure with respect to the second direction. 

4. A method as claimed in claim 1, 2 or 3 wherein the first direction is a 
substantially horizontal direction. 
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5. A method as claimed in any preceding claim wherein the second direction is a 
substantially vertical direction. 



6. A method as claimed in any preceding claim wherein step (e) comprises 
determining p , where 




7. 



A method as claimed in claim 3, or in any of claims 4 to 6 when dependent from 



claim 3, wherein step (d) comprises determining a second-order derivative of the 
pressure with respect to the second direction from the derivative with respect to time of 
the fluid pressure and from the derivative with respect to the first direction of the fluid 
pressure; and using the second-order derivative of the pressure with respect to the 
second direction to determine the first-order derivative of the pressure with respect to 
the second direction. 

8. A method as claimed in claim 7 wherein the second-order derivative of the 
pressure with respect to the second direction is determined using the equation 



claim 3, wherein step (d) comprises deteimining the first-order derivative of the 
pressure in the second direction using the equation: 




9. 



A method as claimed in claim 3, or in any of claim 4 to 6 when dependent from 



°zP\ x q> z q) 1 + 7TA j o^oJ-- 

n 2a 2 



MWo)+o(* 1 ) 
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10. A method as claimed in claim 3, or in any of claim 4 to 6 when dependent from 
claim 3, wherein step (d) comprises determining the first-order derivative of the 
pressure in the second direction using the equation: 



1- 



h 2 (w 2 



\ 



a 



2 22 



J) 



11. A method as claimed in claim 1 wherein step (a) comprises obtaining seismic 
data at two different depths below the surface of the fluid medium. 

12. A method as claimed in claim 1 1 wherein step (a) comprises measuring the fluid 
pressure at two different depths below the surface of the fluid medium. 

13. A method as claimed in claim 12 and comprising determining the first-order 
derivative of the pressure in the second direction using the equation: 

a »fv r 1- g(*o» z o +*/ 2 )- jy(*o.*o -hi 2) h 1 ( \. nfh 4* 



14. A method as claimed in claim 13 wherein step (d) comprises determining d^p 

from a derivative with respect to time of the fluid pressure and from a derivative with 
respect to the ^-direction of the fluid pressure. 

15. A method as claimed in claim 1 wherein step (a) comprises measuring the fluid 
pressure and the particle velocity. 

16. A method as claimed in claim 15 and further comprising the steps of: computing 
the particle velocity from the measured values of the pressure; and comparing the 
results of the computation of the particle velocity with the measured values of the 
particle velocity. 
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17. A method as claimed in claim 15 and further comprising the steps of: computing 
the pressure from the measured values of the particle velocity; and comparing the 
results of the computation of the pressure with the measured values of the pressure. 

18. A method as claimed in claim 15 and comprising the further step of determining 
the depth of the selected one of the seismic receiver and the seismic source from the 
measured values of the pressure and the particle velocity. 

19. A method as claimed in any preceding claim wherein the seismic data is 
obtained in the vicinity of an array of receivers; and wherein the method comprises 
using the results of the processing of the seismic data to reduce the effects of seismic 
energy reflected and/or scattered at the sea surface in the signal emitted by a 
corresponding array of seismic sources. 

20. A method of calibrating a seismic receiver, the method comprising the steps of: 

a) measuring the particle velocity using a seismic source and a seismic receiver 
disposed within a fluid medium; 

b) detennining a derivative with respect to time of the particle velocity at a point in 
the vicinity of the receiver; 

c) detennining a derivative with respect to a first direction of the particle velocity 
at a point in the vicinity of the receiver; 

d) determining a derivative with respect to a second direction of the particle 
velocity at a point in the vicinity of the receiver from the results of steps (b) and (c); and 

e) determining a measure of the quality of the measurement of the particle velocity 
using the derivative with respect to a second direction of the particle velocity. 

21 . A method as claimed in claim 20 wherein step (e) comprises detennining the 
curl of the measured particle velocity. 

22. A method of calibrating a seismic receiver comprising the steps of: 

a) disposing a seismic receiver in a fluid medium; 

b) measuring the fluid pressure and the particle velocity using the seismic receiver; 
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c) computing one of the particle velocity and the fluid pressure from the measured 
value of the other of the particle velocity and the fluid pressure and from the depth of 
the receiver; and 

d) comparing the computed value of the one of the particle velocity and the fluid 
pressure with the measured value of the one of the particle velocity and the fluid 
pressure. 

23. A method as claimed in claim 20 and further comprising measuring the depth of 
the receiver, and using the measured depth of the receiver in step (c). 

24. A method of determining the depth of a seismic receiver, comprising the steps 

of: 

a) disposing a seismic receiver in a fluid medium; 

b) measuring the fluid pressure and the particle velocity using the seismic receiver; 
and 

c) determining the distance between the fluid surface and the receiver from the 
measured values of the particle velocity and the fluid pressure. 
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